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ABSTRACT 

Supermassive black hole binaries (SMBHBs) are expected by the hierarchical 
galaxy formation model in ACDM cosmology. There is some evidence in the 
literature for SMBHBs in AGNs, but there are few observational constraints on 
the evolution of SMBHBs in inactive galaxies and gas-poor mergers. On the 
theoretical front, it is unclear how long is needed for a SMBHB in a typical 
galaxy to coalesce. In this paper we investigate the tidal interaction between 
stars and binary BHs and calculate the tidal disruption rates of stellar objects by 
the BH components of binary. We derive the interaction cross sections between 
SMBHBs and stars from intensive numerical scattering experiments with particle 
number ~ 10 7 and calculate the tidal disruption rates by both single and binary 
BHs for a sample of realistic galaxy models, taking into account the general 
relativistic effect and the loss cone refilling because of two-body interaction. We 
estimate the frequency of tidal flares for different types of galaxies using the 
BH mass function in the literature. We find that because of the three-body 
slingshot effect, the tidal disruption rate in SMBHB system is more than one 
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order of magnitude smaller than that in single SMBH system. The difference is 
more significant in less massive galaxies and does not depend on detailed stellar 
dynamical processes. Our calculations suggest that comparisons of the calculated 
tidal disruption rates for both single and binary BHs and the surveys of X-ray 
or UV flares at galactic centers could tell us whether most SMBHs in nearby 
galaxies are single and whether the SMBHBs formed in gas-poor galaxy mergers 
coalesce rapidly. 

Subject headings: black hole physics — galaxies: kinematics and dynamics — - 
galaxies: nuclei — methods: numerical — X-ray: galaxies 



INTRODUCTION 



In cold dark matter (CDM) cosmology, galaxies form hierarchically and present-day 



galax ies are the products of successive mergers (IKauffmann fc Haehneltll2000l ; ISpringel et al. 



20051 ). Recent observations show that almost all galaxies harbor supermas sive black holes 



(SMBHs) at their centers (IRichstone et al. I Il998l ; iFerrarese fc Fordl 120051 ) and the black 



hole (BH) masses tigh tly correlate with the properties of their ho st galaxies such as the 
mass of stellar bulge (jMagorrian et al" 1998 : Haring fc Rix 2004 ). t he bulge luminosity 
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Gebhardt et al. 



( McLure fc Dunlopl 2002), and the nuclear stellar velocity dispersion ( Ferrarese fc Merritt 



2000l : iTremaine et al.ll2002l ). The correlations between BH mass and 



galaxy properties imply that the growth of SMBHs and the formation and evolution of galax- 
ies are closely linked. The correlation is likely induced by galaxy major mergers (merging 
of galaxies with comparable mass) in which both rapid star formation and gas accretion 
onto SMBHs are triggered and the feedback from the central active galactic nuclei (AGNs) 
regulates the growth of both SMBHs and galaxy bulges. The coevolution scenario can suc- 
cessfully explain not only the correlations between SM BHs and their host galaxies but also 



many of the observed evo lutions of galaxies and AGNs flSpringel et al. 
20061 : ICroton et al.ll2006h . 



20051 : iHopkins et al. 



During galaxy mergers part of the gas originally in gal actic plane is driven to galaxy cen 



ter, t riggering starburst and the accretion of central SMBH ( Gaskelllll985 : Hernquist fc Miho; 



19951 ). If both galaxies in the merging system harbor SMBHs at their centers, a pair of AGNs 
could form which is an intriguing object for both observations and theories. Because of dy- 
namical friction the two SMBHs quickly sink to the common center of the two merging 
galaxies, forming a b ound, compact (~ 10 pc) supermassive black hole binary (SMBHB, 
Begelman et al.lll980l ). which would be difficult to resolve in any but the closest galaxies us- 
ing current telescopes. As the separation of the binary shrinks, dynamical friction becomes 
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less and less efficient and three-body interaction between SMBHB and the stars passing by 
becomes more and more important. At separations ~ 0.1 — 1 pc SMBHB becomes hard in 
the sense that the binary loses energy and angular momentum mainly through three-body 
interaction. The stars passing by the binary, if not tidally dis rupted or swallowed by the 
BHs, would be expelled by the "slingshot effect" (lQuinlanlll996l . Q96 thereafter) with veloc- 
ities comparable to the host galaxy's escape velocity. If the refilling of the reservoir of stars 
on orbits that interact with the SMBHB is inefficient, the evolution of SMBHB slows down 
and the evolution timescale may even exceed the Hubble time, with the SMBHB separation 
stallin g well out side the radii needed for the final gravitational wave radiation and coalescence 
stage (lYul 120021 . Y02 hereafter). This seems inconsistent with that fact that SMBHBs have 
not yet been observed in nearby galaxies. Therefore, many stellar-dynamical mechanisms 
are proposed to boost the hard ening rate of SMBHBs in gas -poor systems, including non- 
equilibrium stellar distr ibution (jMilosavljevic fc Merrittll2003l ). Brownian motion of SMBHB 
( iChatteriee et al.l 120031). and non-spherical (axisymmetric , triaxial, or irregular) stellar dis- 
tribution (Y02; iMerritt &: Poon 112004 iBerczik et al.ll2006l ). However, the conclusions in the 
above scenarios depend on processes which are not well understood, so it is still unclear 
whether SMBHBs could coalesce within a Hubble time using purely stellar-dynamical pro- 
cesses. In gas-rich systems it has been argued that the interaction between binary BHs 
and the gaseous environment could efficiently drive the binary to the gravitational radiation 
dominant domain within a timescale of ~ 10 7 yr (llvanov et al. 



Armitage fc Nataraj an! 120021 ; lEscala et al.ll2005l ; iKazantzidis et al.l 120051 ). 



1999 



Gould k Rix 2000 



Although many physical processes have been proposed in the literature to boost the 
hardening rates of SMBHBs, the abundance of binary BHs in galaxy centers has not yet been 
strongly constrained observationally. Because of their compactness, bound SMBHBs are very 
difficult to resolve directly with current telescopes, a nd only unboun d SMBHBs in a couple of 



gas-rich (wet) galaxy mergers have been identified (IKomossa II2006I ). So far all the evidence 



for bound and hard SMBHBs are indirect and model dependent. The prototype evidence 



19881 : iLiu et al.l Il995t Il997 



emitting BH in SMBHB (IVillata fc Raiteri 



Beeelman et al. 


1980). 


AGNs ( 


Sillanpiia et al. 



Raiteri et al.l 120011) are due to the orbital motion of the jet- 



1999 



and a standard accretion disk (ISillanpaa et al 



dominated accretion flow (ADAF) (ILiu fc Wu 



1988 



Ostorero et al.l 120041 ) , the precession of 



spin axis of the rotating primar y BH jRomero et aL 2000), the interaction between SMBHB 



2002 



Valtaoia et al.l 120001) or an advection 



Liu et al.l 120061 ) . X-shaped features 



in a subclass of radio galaxies have be en attribut ed to the interac tion and alignment of 



SMBHBs and standard accretion disks (ILiu 1 120041 . alternatives see IMerritt fc Ekersl 12002 
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Zier 1120051) . There is also some evidence for the coalescence of SMBHBs in AGNs. iLiu et al. 
( 120031 ) suggested that the double-double radio galaxies and the restarting jet formation in 
some radio galaxies are due to the removal and refilling of inner accretion disk because of 
the interaction with SMBHBs. Although there is much circumstantial evidence for both 
active and coalesced SMBHBs in AGNs, it is still unclear what fraction of SMBHBs would 
coalesce during or before the AGN phase and how many could survive to the later dormant 
or inactive phase. Because present observations cannot give constraints on the physical 
proces ses which have been proposed to boost the hardening rates of SMBHBs, ILiu &: Chen 
( 120071 ) suggested measuring the acceleration of jet precession as a way of determing the 
hardening rates of SMBHBs in AGNs and to place strong constraints on models of SMBHB 
evolution. In inactive galaxies or gas-poor mergers the observational evidence for SMBHBs 
is very rare. One possible way to identify uncoalesced SMBHBs in inactive galaxies is to 
detect hypervel ocity binary stars ejected by SMBHBs with the three-body sling-shot effect 
( ]Lu et al.l 120071 ). However, this method cannot be applied to distant galaxies. Therefore, in 
this paper we suggest a way to determine statistically whether SMBHBs in nearby galaxies 
have coalesced. 

One inevitable impact of a non-rorating SMBH on its stellar environment is that stars 
passing by the BH as close as the tidal radius 

V3 /„ \ / l,T \ -1/3 / M> \ 1/3 



M. 



n — - ~ 10~° P c — 



M 



(lO 8 M 



would be tidally disrupted ([Hills Ill975l ; iRees Ill988l ). where r* and M* are the radius and 
mass of the stars and M. is the BH mass. Part of the stellar debris will be spewed to highly 
eccentric bound orbits a nd later fall back onto the BH, giving rise to an outburst decaying 
within months to ye ars (Rees 19881) . These tidal flares are definitive evidences of SMBHs 



in inactive galaxie s (IKomossa fc Greinerl Il999l ; iKomossa et al.l 120041 ; lHalpern et al.l 12004 



Gezari et al.l 120061 ) . Since the frequency of stellar disruption depends on the surrounding 
stellar distribution (jSyer fc Ulmer I Il999l ; iMagorrian fc Tremaine Ill999l . thereafter MT99), 
the stellar disruption rate can be used as a probe of the inner structure of galactic nucleus. 
For nearby early-type galaxies recent theoretical calculation with single BH without taking 
into account general relativistic (GR) effects gives an averaged stellar d isruption rate of 
~ 10~ 4 yr^ 1 per galaxy if the stellar distribution is spherically symmetric (jWang &: Merritt 



20041 ) or possibly 1 — 2 ord ers of magnitude higher if the stellar distribution is non-spherical 
(axisymmetric or triaxial, iMerritt fc Poon 1 120041 ) . However, if SMBHBs reside in galactic 
centers the stellar distribution would be dramatically cha nged due to t he slin gshot effect, 
and the stellar disruption rate would be very different. Ilvanov et al.l (120051 ) studied the 
interaction between SMBHB and dense galactic cusp (~ 1 pc) for a bound, but non-hard, 
binary at the dynamical friction stage and found a high disruption rate of 10~ 2 — 1 yr _1 . 
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Merritt fc Wang (120051 ) considered the stellar disruption by single BHs immediately after 
SMBHB coalescence. They found that the tidal disruption rate is significantly suppressed 
due to the low density core forming duri ng; the hard stage of SM BHBs. Since the two stages 
investigated by llvanov et al.l (120051 ) and iMerritt fc Wangl (120051 ) are two short phases in the 
evolution of SMBHBs and binary BHs with mass ratios > 10~ 3 spend most of their lifetime 
on the hard stage (Y02), it would be very interesting to calculate the tidal disruption rate 
of stellar objects for hard SMBHBs. Therefore, in this paper we calculate tidal disruption 
rates for hard and stalling SMBHBs. We suggest that the comparison of expected tidal 
disruption rate by single SMBHs or hard SMBHBs and the future survey of X-ray or UV 
flares at nearby galactic centers could tell us whether most of the SMBHs in nearby galaxies 
are single or binary and whether the SMBHBs would coalesce rapidly, as expected by some 
models. 

The stellar disruption rate in a galaxy harboring a SMBHB depends on: (1) the rate 
that stars are fed from large-scale to the vicinity of SMBHB and (2) the probability that 
a star moving in the vicinity of the SMBHB passes close enough to one of the BHs to 
be tidally disrupted. Although there has been much work on the former part, inte nded 
to study the evolution timescale of SMBHB (Q96; Y02; iMilosavlievic fc Merrittll2003h . the 
latter is largely unaddressed in the literature. Because of the chaotic nature of three-body 
interaction, numerical scattering experiments are needed to give the probability of tidal 
disruption. Previous numerical experiments on star- SMBHB inter actions mainly focus on 
energy and angular momentum exchanges (Q96; ISesana et al.ll2006l ) and the limited particle 
numbers are insufficient to tackle the rare events of very close encounters such as tidal 
disruptions. Therefore, we begin with intensive numerical scattering experiments (normally 
~ 10 7 particles in each run) to calculate the tidal disruption cross sections for hard SMBHBs. 
Then we apply our results to a sample of nearby galaxies and estimate the tidal disruption 
rates of stellar objects both by single and binary BHs in different types of galaxies. 

This paper is organized as follows. In § [2] we briefly review some basic properties of 
hard SMBHBs and show how to calculate stellar disruption rate. We describe our numerical 
scattering experiments in § [3] and present the results in § HI In § [5] we calculate stelar 
disruption rates with realistic galaxy models and discuss the differences between the tidal 
disruption rates for single and binary BHs. We also estimate the density of tidal disruption 
events in local universe and make predictions for future surveys of UV/X-ray flares in § [6j 
In § [7] we discuss the observational signatures of SMBHBs in inactive galaxies and the 
implications for the evolution of binary BHs. Finally we give our conclusions in § [SI 
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2. PROPERTIES OF HARD SMBHB SYSTEMS 

2.1. Loss Cone Filling Rate 

A SMBHB with masses Mi and M 2 {Mi > M 2 ) becomes hard when the semimajor axis 
a decreases to 

a ~ ah = 4^2 > ( 2 ) 

(Q96) where is the hardening radius, cr* is the one-dimensional stellar velocity disper- 
sion of background stars, and \i = M\M 2 lM\ 2 is the reduced mass with M 12 = Mi + M 2 . 
Observatio ns of nearby galaxies show that BH mass tightly correlates with bulge velocity 



dispersion flFerrarese fc Merrittl l2000l ; iGebhardt et al. 1120001 ). If the total BH mass M 12 of 
SMBHB also follows the empirical M, - a* relation: M u = lO 813 Af (a*/2OO km s" 1 ) 402 



(ITremaine et al.ll2002l ). equation (j2J) reduces to 



a fe ~3.13 P c - q X9 M 8 1/2 , (3) 



where M§ = Mi 2 /10 8 M Q and q = M 2 /Mi is the mass ratio of the binary. In this paper we 
only consider q ^ 0.01 because the dynamical friction timescale for SMBHBs with q 0.01 
is very long (Y02). 

A hard SMBHB loses energy and becomes more bound through three-body interaction 



f lHut II1983I ). As a result stars intruding a sphere of radius ~ a about the mass center of the 
binary (if not tidally disrupted or swallowed by the BHs) will be eventually expelled with 
an average energy gain 

GMi 2 M*Aa , , . 

= KGfiM./a , (4) 

where A is a dimensionless factor about 1.6 (Q96; Y02). In this paper, unless noted other- 
wise, we assume that galaxies are spherical. In a spherical potential there is just one family 
of orbits (centropilic loops) and any such orbit can be labelled by its four integrals of motion, 
namely the specific angular momentum J and specific binding energy S. We define the loss 
cone to consist of all orbits with pericenters less than a, which results in a conical region 
with the boundary given by J\ c (£, a) = 2a 2 [GMi 2 /a — £]. In a non-spherical galaxy angular 
momentum is no longer conserved and there can be more than one orbit family. Nevertheless, 
instead of a "loss cone" one can still define a "loss zone" inside which stars will be delivered 
to the SMBHB. 



Mechanisms such as two-body relaxation tend to refill the loss cone so the decay of 
SMBHB continues. The stars refilled into loss cone are originally far from the SMBHB so 
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that 



and 



8 < GM 12 /a 



(5) 



2GM 12 a . 



(6) 



If the characteristic change of J during one orbital period of a star is much smaller than Jj c , 
the loss cone refilling behaves like a di ffusion process in which the loss cone filling rate is not 
sensitive to Jf c or a (diffusive regime, ICohn fc Kulsrudlll978l ). Otherwise, the stars entering 
the loss cone have a uniform distribution with respect to J 2 so that the loss cone filling rate 
is proportional Jf c or a (pinhole regime). 

In inactive galaxies gas dynamics is probably not important, so a hard binary loses 
energy mainly through interacting with loss cone stars. The evolution timescale th at this 
hard stage depends on the loss cone filling rate F lc (a), the number of stars fed into a sphere 
of radius a about the mass center of the binary per unit time. From equation (j3J) and 
d(GM 12 /2a)/dt = F lc (a)AE we have 



t h = \a/a\ = M 12 /[2KF lc (a)m*] 



(7) 



(see also Y02), which implies that about an amount of F lc (a)th ~ M\ 2 /2K of sta rs should 



be co nsumed to reduce a by a single e-fold, in agreement with recent simulations (IMerritt 
2006h . 



2.2. Tidal Disruption Cross Section 

A fraction of the stars fed into the loss cone of SMBHB would be scattered to the 
vicinity of a BH and get tidally disrupted. The ratio of the tidal radius of the primary BH 
and the separation of the binary is about 



^ ~ 3.71 x 10- 5 M- 2/3 I ^ 



M* V 1/3 (l + <?) 5/3 / o 



a h \r Q J \M Q J q_ x V200 km s" 



3.20 x 10~°M 8 ~ VU ( — | 



re J \M Q J q-i 

(from eq. [I] and [3]), which is small for typical hard SMBHBs This implies that most stars 
in the loss cone would be expelled instead of being tidally disrupted. 



To further quantify the probability of tidal disruption, we introduce Sj(r) (i = 1,2), 
the cross section for stars to be scattered to a distance less than r from the ith BH. The 
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probability that a loss cone star is disrupted by the ith BH estimates £j(r t j)/Ej(a), where 
r ti is the tidal radius of the ith BH. Then the stellar disruption rate is related to the loss 
cone filling rate and the tidal disruption cross section as 



For a single BH the geometrical cross section can be written out analytically as 

™-*( 1 + g &)' (10) 

which increases with the gravitational potential GM,/r and decreases with stars' kinetic 
energy v 2 /2 at infinity. For bound stars equation (fT0|) is still valid as long as the stars fall to 
the BH on near-radial orbits, and in the bound case Vq is the velocity at apocenter. When 
gravitational focusing is important (GM,/r 3> v 2 /2) the cross section scales as r. In SMBHB 
systems the presence of a companion BH tends to increase the term GM,/r by deepening 
the gravitational potential and also increase the term Vq/2 by inducing orbital motion which 
raises the relative velocity between each binary member and ambient stars. Therefore in 
the rest frame of the primary member, the equivalent increase in GMi/r is GM 2 /r and the 
equivalent increase in Vq/2 is about q 2 GMi2/a. So the cross section of the primary BH is 
approximately 

~ , s •> ( GM 12 /r \ , , 

where 77 is a correction factor of order unity. We will see below that for different qs the best 
fit is obtained when 77 = 1/(1 + q). For the secondary, the roles of Mi and M 2 exchange, 
so the cross section can be obtained by replacing q in equation (lllj) with 1/q. Quinlan has 
separately explained the behavior of the cross sections at r/a ^> q and r/a <C q from another 
way of understanding (Q96), while our explanation and empirical formula are general for any 
r/a. 

Because of the effect of GR equation (iTTjl is no longer valid when r ~ rs, where rg = 
2GM,/c 2 is the Schwarzschild radius. Since the ratio of r$ and r t is 



GR effect is important when M, > 3.6 x 10 6 M© (r t < 10rs for solar type stars). In practice, 
one can use the pseudo-Newtonian potential 



- 9- 



( jPaczyriski fc Witta Ill980l ) to simulate the GR effect so equation ffTTl) becomes 

GM 12 /(r-r sl ) 



Sir 



7rr 



f + 



w 2 /2 + m 2 GM 12 /a 



(14) 



(notice that a 3> r Si ), equivalent to multiplying the non-relativistic cross section by a correc- 
tion factor r/(r — rg). Equation (1141) diverges when r — > rs, but since a star plunging into the 
sphere of marginally bound radius 2rs about a BH eventually falls into to the event horizon, 
the cross section at [rs, 2rs] should be constant, equal to the cross section at r = 2rs. 

Both equation ffTTj) and (fT4"l) are evaluated at the first time that a star passes the 
binary. In some cases a star could be scatted onto tempor ally bound orbit and encounters 
with the binary many times before expelled or disrupted (IHills Ill983l ; Q96). During such 
multi-encounters the probability of tidal disruption is expected to be enhanced. Since it is 
very difficult to derive analytical cross sections for these complicated encounters, scattering 
experiments are needed. Stars on near-radial orbits have Vq <C GM^/a, then equations (IE]), 
( TTTj) . and (fT4"j) suggest that the probability of tidally disruption scales as r ti /a, which is of 
order 10 -4 — 10 -6 . Therefore, to get statistically meaningful tidal disruption cross sections 
a large number of particles should be used in the scattering experiments. 



3. METHOD 
3.1. Scattering Experiments 

In this section we describe our numerical scattering experiments of the restricted three 



body problem. The method adopted here is similar to those in Q96 and lSesana et all (120061 ) 
but the number of test particles in our experiments is much larger, typically ~ 10 7 , to give 
statistically meaningful cross sections at 10 -6 < r/a < 10 -4 . 

In the scattering experiments, instead of using S, the specific binding energy of a particle 
in the combined potential of stars and BHs, we use a more convenient parameter E = 
GM\ 2 jr — v 2 /2, the specific binding energy about the BH binary, where r is the radius 
about the mass center and v is the velocity of particle. We denote the initial specific binding 
energy with Eq and initial velocity with vq. In our problem the relevant energy range is 
-qGM 12 /Sa < £ < GM 12 /a (from eq. [2] and eq. p]). 

In each scattering experiment, the origin is chosen at the mass center of the binary. 
In the case —qGMi 2 /8a < Eq < (unbound case), initially particles come form infinity 
r = (x,y,z) = (b, 0, oo) with asymptotic velocity v = (v x ,v y ,v z ) = (0,0, —v ), where b 
is the impact parameter, the minimum separation between particle and mass center if the 
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particle feels no gravitational field. The fiducial value for Vq is 10~ 3 ^GMi 2 /a, to reproduce 
Fig. 3 in Q96. Given t> , b 2 is uniformly sampled in the range [0,6^ ax ], where b max = 
y2Jic(a)/vo- Particles with b > 6 max hardly reach r < a therefore contribute little to the 
tidal disruption cross section (see § I4.3| . In the other case ^ E <C GM^/a (bound case), 
initially we put particles at r = (x, y,z) = (0,0, r ) with an isotropic velocity distribution 
(t> = a/2(GMi2/Vo — E )). r satisfies a « r < GM^/Eq because in realistic SMBHB 
systems most stars enter the loss cone on near-radial orbits. 

Another six parameters (four if the binary is circular) are set to fix the initial conditions 
of the binary: (1) the mass ratio q; (2) the eccentricity of binary orbit e; (3) the inclination 
of orbital plane; (4) the argument of pericenter; (5) the longitude of ascending node; (6) the 
initial binary phase. In each set of experiments with fixed q and e, the cosine of the orbital 
inclination angle is evenly sampled in [—1, 1] and the other three angular parameters are 
uniformly distributed in [0, 27r] , resulting in an isotropically filling of the loss cone. 

Given the above initial conditions, the motion of a massless particle is governed by the 
following coupled first-order differential equations: 



. = ^m^iA (us) 

^— ' r — Yi \ 6 

i=i 

where is the position of the ith BH. In each scattering experiment we first move a particle 
from its initial position to = a(10 8 g) 1 / 4 along a Keplerian orbit about a point mass Myi 
at the origin. At the quadrupole force from the binary is eight orders of magnitude 
smaller than GM^o /a 2 . Then we integrate the particle's orbit with the subroutine dopri8 



( iHairer et all 119871 ). an explicit Runge-Kutta method of order (7)8. We set the threshold 
of fractional error per step in r and v to be 10 -8 . Raising this threshold up to 10~ 6 does 
not significantly change our results. The integration stops if one of the following conditions 
is satisfied (1) the particle leaves the sphere of radius with negative binding energy; (2) 
the physical integration timescale exceeds 10 10 yr; (3) the integration timestep reaches 10 6 . 
A small fraction (< 0.1 percent) of the particles are scattered to wide, bound orbits and 
survive many revolutions so conditions (2) and (3) are adopted to save computational time. 
We have tested our code by reproducing Fig. 3 in Q96 and found very good agreement. 

We have also used the pseudo- Newtonian potential (eq. [T3]) in our numerical experi- 
ments to investigate the effect of GR. In this case equation (fT5|) becomes: 

M t (r-r,) 

2 ' 



— 
s-^i r _ r 



i=l 



\r-Ti\- r Si ) 
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and we stop the integration once a particle reaches 1.01 Schwarzschild radius about either 
of the BHs. In the GR experiment the ratio of rsi and ah, 

^ ~ 3 .0xl0- 5 M 8 1/2 i±^, (16) 
a h 8 q-i 1 ^ 

should be set in priori. For the illustrative purpose we always set M§ = 1. We will discuss 
the effect of changing M 8 on tidal disruption cross sections in § 14.51 



3.2. Derivation of the Cross Sections 

After each scattering experiment we record the minimum separation between the par- 
ticle and each BH. At the end of all experiments we count the number iVj(r) of parti- 
cles whose minimum separations from the ith BH are less than r. Then the normalized 
multi-encounter cross section (particle are allowed to encounter the binary as many times 
as they could until they are expelled) is calculated with £j(r)/Ej(a) = A^(r)/A^(a). The 
Poissonian error in the counts is yiV^r), so the error for the normalized cross section 
is Ni(r)/Ni(a)y/l/Ni(r) + 1/Ni(a) and the corresponding fractional error in statistics is 
<Wr) = y/l/Ni(r) + l/Ni(a). 

Following Q96, we also record the minimum separation between a particle and each 
BH during their first encounter and derive the single-encounter cross sections. In § 14.11 we 
compare the single-encounter cross sections and those from the analytical approximations 
(eq. P3] and [H]). 



4. RESULTS FROM THE SCATTERING EXPERIMENTS 

4.1. Non-relativistic Cross Sections 

First we consider the unbound case. We set vq = 10~ 3 ^/ GM 12 /a (after Fig. 3 in 
Q96), q = 0.01, and use the Newtonian potential. The normalized mullti- and single- 
encounter cross sections from N ~ 10 7 particles are presented as solid lines in the top 
panel of Figured! The cross sections are plotted as a function of r/a so that they can be 
easily scaled. The perturbations at the lower ends are due to statistical fluctuation. In 
the top panel of Figure [1] we also plot the empirical single-encounter cross sections (see 
eq. |llj ) as dashed lines. We estimate the fractional error of the empirical formula with 
Oap P = |£i(r)/£j(a) - S i (r)/S i (a)|/[S i (r)/S i (a)], where Sj(r)/Sj(a) is from scattering ex- 
periments and £i(r) from equation ( fill . cr app and 0" sta t for the primary and secondary BHs 
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are shown in the middle and bottom panels of Figure [TJ For both BHs we found that at 
r <C a (T a pp is almost always below <7 sta t, indicating that the empirical cross sections agrees 
very well with the numerical results. At r > qa, although cr app is larger than cr sta t, the 
difference does not exceed 10%. 

For q = 0.1 and q — 1, the numerical (solid lines) and empirical approximations (dashed 
lines) cross sections are presented in the top left and top right panels of Figure [2j cr stat 
(solid lines) and <j app (dotted lines) for the single-encounters are presented in the middle and 
bottom panels. Although the approximations become worse when q is large, the difference 
between approximation error and statistic fluctuation is always below 10%. 

Figure [1] and Figure [2] show that at r < q 2 a both the multi- and single-encounter cross 
sections scales linearly as r/a, implying that gravitational focusing is important for very 
close encounters. The cross sections at r/a < 10~ 6 could be obtained by using this scaling. 
It is also clear that the normalized multi-encounter cross sections are greater than the single- 
encounter ones as is predicted in § 12.21 However, we find that in our experiments although 
a particle could encounter with the binary many times, no particle more than once enters 
the sphere of tidal radius (for M 8 = 1) about either BH. This is because the probability 
for a particle to enter the sphere of tidal radius about the ith BH n times roughly scales 
as [£j(rtj)/Ej(a)] ra or (r t i/a) n , which becomes extremely low for r t i <C a and n > 1. For 
the similar reason, the event that a particle successively approaches the tidal radii of the 
primary and secondary BHs is also extremely rare. 

4.2. Dependence on Particle's Binding Energy 

To investigate the effect of binding energy, we carried out two test experiments for 
Eq = —0.01 and Eq = 0.01 (in unit of GM^/a), each with ~ 10 6 particles and q = 0.01. In 
the unbound case the asymptotic velocity of the intruding particles is v ~ 0.14-^/ GM^/a, 
slightly greater than the limit 0.85 -^/GM^q/ a (from the fitting formula eq. 17 in Q96) above 
which the binary becomes soft. In the bound case, initially the particles are at r = 50a 
with an isotropic velocity of Vq ~ 0.14^/GM 12 /a. We do not consider Eq = 0.1 because such 
initial condition results in unrealistically bound particles with tq ^ 10a. 

Results from these test experiments are presented in Figure [3] as solid lines. Cross 
sections for the fiducial value Eq ~ — 10 _6 GM 12 /a (from § !4.ip are also plotted for reference. 
Although the physical cross sections increase prominently with Eq which is expected by 
equation [Til Figure [3] shows that in both bound and unbound cases the normalized cross 
sections for both multi- and single-encounters seem not varying significantly with Eq. This 
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Fig. 1. — Top: Normalized multi- and single-encounter cross sections from 10 7 particles for 
q = 0.01. From top to bottom, the first two solid lines, respectively, refer to the multi- 
and single-encounter cross sections of the primary BH, and the third and fourth ones are 
of the secondary. Empirical cross sections (dashed lines) and positions of tidal radii and 
Schwarzschild radii for M\2 = 1O 8 M (short vertical lines ) are also plotted. Middle: Frac- 
tional errors of statistical fluctuation (solid) and empirical cross sections (dotted), for the 
primary BH. Bottom: The same as the middle panel but for the secondary BH. 




Fig. 2. — Close-encounter cross sections and fractional errors for q — 0.1 (the left panels, 
10 7 particles) and q = 1 (the right panels, 10 6 particles). Lines have the same meanings as 
those in Figure [U 
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Fig. 3. — Normalized multi- and single-encounter cross sections from ~ 10 6 particles (solid 
lines, q = 0.01) for E = 0.01 (top) and E = -0.01 (bottom). E is in unit of GM 12 /a. 
Dotted lines are from the fiducial experiments (§ 14.11) and short vertical lines have the same 
meanings as those in previous figures. 
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is because when \Eq\ <C GM^/a, the velocity of the particles passing r < a from the binary 
is not sensitive to Eq but to the depth of the gravitational potential at r. Therefore, different 
particles passing by the same BH binary feel the similar strength of gravitational focusing 
and have similar interaction cross sections with the BH components of binary. Since the 
normalized cross sections do not vary significantly with Eq as long as \Eq\ <C GM^/a, it 
is reasonable to apply the fiducial cross sections to various binaries with a wide range of 
semimajor axis. For example, according to equation ([2]), varying Eq from — 10~ 6 GM 12 /a to 
— 10~ 2 GM 12 /a corresponds to increasing a from ~ lQ~^qZ\a h to ~ O.lglja^. 



4.3. Dependence on Particle's Initial Angular Momentum 

In both experiments for Eq ^ and Eq < initially the particle follow a uniform 
distribution in J 2 and the loss cone filling rate is isotropic, i.e., the cross sections obtained in 
§ 14.11 and § 14.21 have been equally averaged over the whole loss cone and over all directions. 
Such averaging may be oversimplified in realistic SMBHB systems. 

To investigate the dependence of the cross sections on the amplitude of initial angular 
momentum, in FigureHJwe plot the differential tidal disruption cross sections dEj(r^) / 'dJ 2 /E,(a) 
with respect to J 2 . The triangles and squares, respectively, refer to the primary and sec- 
ondary BHs. dHi(r t i) I 'dJ 2 /Ej(a) is obtained as follows. For each bin of A J 2 we count 
the number AA^(r^) of particles whose initial angular momenta fall in this bin and min- 
imum separations from the zth BH are less than r t i- Here r t i is from equation (|Sj) with 
M 8 = 1 and we only use the data for multi-encounters. Then dTs/dJ 2 is calculated with 
ANi(r ti )/AJ 2 /Ni(a) so that the integration over J 2 results in the normalized cross section 
of tidal disruption. The Poissonian errors in the counts have been indicated by the error bars 
in Figure HI Note that the differential cross section [or AiV* (?■*«) alone] derived in this way is 
proportional to the probability of tidal disrupting a particle with initial angular momentum 
J. Generally speaking, both differential cross sections for the primary and secondary BHs 
are flat at J 2 < J 2 C ~ 2GM 12 a. The differential cross sections start to cut off above Jf c and 
particles with initial angular momentum greater than 2Jj 2 contribute little to the cross sec- 
tions. When q <C 1 the differential cross section for primary BH rises steeply with decreasing 
J 2 at J <C Ji c but even in this case the contribution to the total cross section from particles 
with larger J is still significant. 

The cross sections also depend on the direction of the initial angular momentum vector 
J. We define 9 as the relative angle between J and binary's orbital angular momentum. 
Figure [5] shows the differential tidal disruption cross sections with respect to cos 6. The 
differential cross sections are obtained following the same described above. When q <C 1 the 
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Fig. 4. — Differential tidal disruption cross sections with respect to the square of initial 
angular momentum (in unit of GM 12 a) for q = 0.01 {top), q = 0.1 (middle), and q = 
1 (bottom). Triangles and squares, respectively, refer to the differential cross sections of 
primary and secondary BHs. Error bars indicate the Poissonian errors in the counts of tidal 
disruption events. When differential cross section is low, large perturbation occurs due to 
the noise in the counts. 
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particles on corotating orbits (cos 6 > 0) are more likely to be disrupted, while for secondary 
BHs the dependence of the cross sections on 9 is weak. Unless the stellar distribution is 
extremely flattened or the star cluster around SMBHB is significantly rotating the total 
cross section would not differ much from those presented in Figures [T] and [2J 

4.4. Effect of Eccentric Binary Orbit 

To see clearly the effect of the eccentricity of binary's orbit, we set an extreme eccen- 
tricity e = 0.5 and carried out two test experiments for q = 0.01 and 0.1, each with 10 6 
particles. The results (solid lines) are presented in Figure El Cross sections for e = (dotted 
lines, from § 14.11) are also plotted for reference. Increasing eccentricity only slightly increase 
the normalized multi-encounter cross sections for both BH components. 

4.5. Effect of GR 

To study the effect of GR we set the gravitational potential to be pseudo-Newtonian 
with equation [13] and Mg = 1 and then repeated the scattering experiments. For each q the 
result from 10 7 particles is shown in Figure [7] as solid lines. Dotted lines are results from the 
non-relativistic experiments. Below the marginally bound radius 2rg both the cross sections 
for multi- and single-encounters become constant, representing the relativistic effect. 

The dashed lines are empirical cross section in the GR case. For the single-encounters 
we have used equation (TTJJ and set the cross sections constant between r = rs and 2rs- The 
resulting cross sections agree well with the numerical ones. The approximate cross sections 
to the multi-encounters are obtained by multiplying the non-relativistic cross sections by 
the correction factor r/(r — rs) and keeping the cross sections constant below 2r$. These 
resulting cross sections well agree with the numerical ones as long as rs <C q 2 a, i.e., in the 
gravitational focusing dominated regime. Due to GR effect the cross sections at r = rs are 
four times greater than the non-relativistic ones. 

For other values of M 8 the location of Schwarzschild radius moves along the r/a axis. 
Instead of repeating the whole simulation we use non-GR cross sections to calculate the 
approximate relativistic cross sections, following the procedure described above. 
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Fig. 5. — Differential tidal disruption cross sections with respect to cos 9 for q = 0.01 (top), 
q — 0.1 (middle), and q = 1 (bottom). Triangles and squares, respectively, refer to the 
primary and secondary BHs. 
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Fig. 6. — Normalized multi- and single-encounter cross sections from 10 6 particles (solid 
lines) for elliptical binaries with e = 0.5, q = 0.01 (top) or q = 0.1 (bottom). Dotted lines are 
results from the circular binary case. Short vertical lines have the same meanings as those 
in Fig. [TJ 
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Fig. 7. — Relativistic cross sections from 10 7 particles (solid lines) for q = 0.01 (top left), 
q — 0.1 (bottom left), and q = 1 (bottom right). The empirical cross sections in GR case 
(dashed lines) and the cross sections from non-GR experiments (dotted lines) are also plotted. 
Short vertical lines have the same meanings as those in Fig. [TJ 
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STELLAR DISRUPTION RATES IN NEARBY GALAXIES 



5.1. Galaxy Sample 



In this section we estimate the tidal disruption rates (eq. ED of stellar objects for a 
sample of nearby galaxies. The largest uncertainty in our calculation comes from the loss 
cone refilling rate F lc (a); for simplicity, and to ease comparison with earlier work, we assume 
that two-body relaxation is the only process contributing to F lc (a). 

Our galaxy sample consists of 51 nearby elliptical galaxie s which are listed in Table [TJ 
For each galaxy the BH mass is obtained with lHaring fc Rix r s M. — Mb u i ge relation, where 
the bulge mass Mh„i g e is from the scaling relation between stellar mass and luminosity 
( iMagorrian et al.lll998l ). The assumed mass to light ratios and the resulting BH masses 
are given by columns 2 and 3 in Table [TJ 



Table 1. Galaxy sample 





T v 


log Af. 


log F 1 ? , 

lJ Single 


l°g a stall 


log *evol 




log N? 


logiVf 


logiVf 


log iV^ 


logATf 


logiV| 


Name 


\ XL/ / \ J J 


(M Q ) 


(yr- 1 ) 


(pc) 


(yr) 


(yr" 1 ) 


(yr" 1 ) 


(yr" 1 ) 


(yr" 1 ) 


(yr" 1 ) 


(yr- 1 ) 


(yr" 1 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(H) 


(12) 


(13) 


A2052 


7.41 


9.17 


-5.89 


1.76 


13.55 


-4.88 


-9.75 


-12.54 


-10.33 


-11.87 


-11.11 


-11.41 


NGC 1023 


4.88 


7.84 


-4.13 


0.25 


9.98 


-2.65 


-6.47 


-9.71 


-7.03 


-8.54 


-7.59 


-7.61 


NGC 1172 


5.39 


8.16 


-4.21 


1.09 


10.50 


-2.85 


-7.38 


-10.20 


-7.97 


-9.49 


-8.76 


-9.06 


NGC 1316 


7.60 


9.25 


-4.64 


0.89 


12.46 


-3.71 


-7.72 


-10.47 


-8.25 


-9.80 


-9.04 


-8.94 


NGC 1399 


6.33 


8.67 


-5.28 


0.87 


12.73 


-4.57 


-8.73 


-11.60 


-9.27 


-10.83 


-10.10 


-10.10 


NGC 1400 


5.69 


8.33 


-4.79 


0.97 


11.49 


-3.67 


-8.05 


-11.33 


-8.59 


-10.17 


-9.57 


-9.57 


NGC 1426 


5.04 


7.95 


-4.44 


0.95 


10.38 


-2.94 


-7.45 


-10.29 


-7.97 


-9.51 


-8.78 


-8.78 


NGC 1600 


7.44 


9.19 


-5.80 


1.33 


13.81 


-5.13 


-9.57 


-12.49 


-10.13 


-11.70 


-10.97 


-10.91 


NGC 1700 


6.25 


8.63 


-4.37 


1.18 


10.84 


-2.71 


-7.22 


-10.06 


-7.74 


-9.24 


-8.55 


-8.55 


NGC 221 


2.71 


5.96 


-3.68 


-3.34 


8.52 


-3.16 


-4.96 


-6.48 


-4.84 


-6.07 


-5.37 


-5.37 


NGC 224 


4.62 


7.67 


-4.54 


-0.47 


10.91 


-3.75 


-7.00 


-9.81 


-7.49 


-8.98 


-8.04 


-8.10 


NGC 2636 


3.93 


7.15 


-5.05 


0.88 


9.95 


-3.30 


-7.99 


-11.30 


-8.53 


-10.10 


-9.49 


-9.49 


NGC 2832 


7.76 


9.32 


-5.53 


1.57 


13.57 


-4.75 


-9.41 


-12.69 


-9.94 


-11.44 


-10.95 


-10.95 


NGC 2841 


4.66 


7.69 


-4.59 


0.74 


10.21 


-3.02 


-7.36 


-10.68 


-7.92 


-9.46 


-8.92 


-8.92 


NGC 3115 


5.39 


8.16 


-4.03 


0.56 


10.19 


-2.53 


-6.56 


-9.59 


-7.12 


-8.62 


-7.87 


-7.80 


NGC 3377 


4.53 


7.60 


-3.85 


-0.40 


10.04 


-2.94 


-6.28 


-9.07 


-6.77 


-8.26 


-7.37 


-7.36 


NGC 3379 


5.21 


8.05 


-4.93 


0.63 


11.48 


-3.93 


-8.07 


-10.89 


-8.62 


-10.13 


-9.31 


-9.31 


NGC 3599 


4.55 


7.61 


-4.61 


0.92 


10.16 


-3.05 


-7.63 


-10.41 


-8.18 


-9.70 


-9.04 


-9.26 


NGC 3605 


4.13 


7.31 


-4.78 


0.44 


10.00 


-3.19 


-7.36 


-10.26 


-7.94 


-9.42 


-8.79 


-8.87 


NGC 3608 


5.48 


8.21 


-4.79 


0.94 


11.12 


-3.41 


-7.80 


-11.07 


-8.35 


-9.91 


-9.31 


-9.31 


NGC 4168 


6.39 


8.70 


-5.56 


1.20 


12.77 


-4.58 


-9.09 


-11.93 


-9.61 


-11.11 


-10.45 


-10.39 


NGC 4239 


3.50 


6.78 


-5.24 


-0.02 


10.02 


-3.74 


-7.67 


-10.47 


-8.19 


-9.70 


-8.90 


-8.90 


NGC 4365 


6.71 


8.86 


-5.24 


0.96 


12.82 


-4.46 


-8.65 


-11.50 


-9.21 


-10.69 


-10.06 


-10.14 


NGC 4387 


3.95 


7.17 


-4.83 


0.06 


10.02 


-3.36 


-7.26 


-10.09 


-7.78 


-9.25 


-8.40 


-8.49 


NGC 4434 


4.01 


7.22 


-4.48 


-0.37 


9.98 


-3.27 


-6.74 


-9.42 


-7.25 


-8.74 


-7.86 


-7.85 


NGC 4458 


4.01 


7.22 


-4.31 


-0.39 


9.99 


-3.28 


-6.73 


-9.43 


-7.25 


-8.73 


-7.86 


-7.84 


NGC 4464 


3.54 


6.82 


-3.99 


-2.38 


9.78 


-3.56 


-5.73 


-7.62 


-5.79 


-7.13 


-6.40 


-6.40 


NGC 4467 


2.92 


6.20 


-4.52 


-2.89 


9.59 


-3.99 


-6.01 


-7.73 


-5.96 


-7.27 


-6.55 


-6.55 


NGC 4472 


7.29 


9.12 


-5.30 


0.83 


13.29 


-4.67 


-8.66 


-11.43 


-9.20 


-10.72 


-9.95 


-9.89 


NGC 4478 


4.49 


7.58 


-4.78 


0.80 


10.27 


-3.20 


-7.66 


-10.55 


-8.20 


-9.76 


-9.03 


-8.97 


NGC 4486 


7.05 


9.02 


-5.57 


0.83 


13.40 


-4.89 


-7.27 


-10.08 


-7.76 


-9.25 


-8.34 


-8.40 


NGC 4486b 


3.18 


6.48 


-5.11 


-0.89 


10.00 


-4.03 


-8.91 


-11.65 


-9.45 


-10.98 


-10.24 


-10.14 


NGC 4551 


4.10 


7.28 


-4.70 


0.01 


10.00 


-3.23 


-7.03 


-10.27 


-7.56 


-9.09 


-8.18 


-8.19 



to 
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5.2. Loss Cone Filling Rates 

For ea ch galaxy we buil t simp le spherical models using the same Nuker law parameters 



as those in IWang fc Merrittl (120041 ) and separately calculated the loss cone filling rates due 
to two-body interaction in single and binary BH cases (see MT99 and Y02 for detailed 
description of the calculation) . 

For single BHs we did two sets of calculations to investigate the effect GR on loss cone 
filling rate F^^ e . In the non-GR case we calculated J\ c with r t from equation ([!]) while in 
the GR case we substituted r t with 2r s if r t < 2r s . In both cases the resulting loss cone 
filling rates are similar, within a factor of 1.4 even when 2rs ^> r t . This is because the loss 
cones of the most massive BHs are diffusive so that i^gie is insensitive to the size of loss 
cone. In Tabled] we only give i^^e for the non-GR case. The mean of log-F^f le is —4.82 
for all the galaxies or —4.44 for galaxies with M, < 1O 8 M , which is within a factor of 2 



of the rates obtained by lWang fc Merrittl (120041 ). the differences being completely accounted 



for by our different assumed BH masses. 

For binary BHs GR effect is not important to loss cone filling rate because we are 
interested in the case a 3> r$- For each galaxy we have calculated the loss cone filling rates 
for binary BHs with a set of a and derived the evolution timescale t cvo \(a) with 

1 _ 1 | 1 ^ 

£ e voi(«0 t h (a) t p (a) 



(Y02), where t gr is the evolution timescale due to gravitational wave radiation (IPeters Ill964l ). 
Then the binary's stalling radius a sta ii is determined by locating the peak of the t evo \(a) curve. 
If tevoi(ostaii) > 10 10 yr, which is the case for most of our model galaxies, we increase a sta ii 
until (1) tevoi(ostaii) < 10 10 yr or (2) a sta n > 0.1 arcsec, whichever comes first. The latter 
criteria is motivated by that SMBHBs in nearby galaxies, if wider than 0.1 arcsec, would 
have been observed. Notice that a sta ii, £evoi(<istaii)> an d i 7llc (a s taii) given by either criteria (1) 
or (2) are independent of q. The results for q = 0.1, the typical mass ratio in the literature 
for the present-day SMBHBs, are given in Table (U while the results for q = 0.01 and 1 are 
essentially identical. Although the loss cone of SMBHB in Eq — J 2 space could be more 
than four orders of magnitude wider than that of a single BH (eq.s [6] and [8]), the loss cone 
filling rate for binary BH increases only by a factor of 10, which implies that the loss cone 
of SMBHB is also in the diffusive regime. 



Table 1 — Continued 





TV 


lno- A/7" 


lno- F lc 

10 S ^single 


10 5 tt stall 


l°g ^cvol 


lno F lc 
10 S ^binary 


lno- !V d 
log iVj 


log iV 2 


lno- N d 
log 


log iV 2 


lno- N d 
log iV-j^ 


lno- N d 
log JV 2 


Name 


1 \fl IT \ 






\V C ) 


ly r ) 


(y r ) 


vy r ) 


(y r ) 


ly r ) 


(y r ) 


(y r ) 


ly r ) 


i\\ 
('-) 


( 2 ) 


(■>) 




(.5) 


{°) 




l»J 


/'o^ 


(1U) 




(12) 


(IS) 


NGC 4552 


5.66 


8.32 


-4.81 


0.82 


11.79 


-3.98 


-8.19 


-11.03 


-8.77 


-10.32 


-9.58 


-9.88 


NGC 4564 


4.72 


7.73 


-4.24 


0.23 


10.00 


-2.77 


-6.61 


-9.86 


-7.16 


-8.66 


-7.79 


-7.88 


NGC 4570 


4.80 


7.79 


-4.11 


-0.20 


9.97 


-2.69 


-6.16 


-8.93 


-6.67 


-8.14 


-7.25 


-7.23 


NGC 4621 


5.88 


8.43 


-4.11 


0.82 


10.53 


-2.60 


-6.78 


-9.65 


-7.35 


-8.83 


-8.14 


-8.14 


NGC 4636 


6.28 


8.65 


-5.52 


0.80 


12.95 


-4.81 


-8.93 


-11.85 


-9.49 


-11.00 


-10.14 


-10.14 


NGC 4649 


6.80 


8.90 


-5.37 


0.82 


13.09 


-4.70 


-8.75 


-11.76 


-9.30 


-10.79 


-10.04 


-9.97 


NGC 4697 


5.64 


8.30 


-4.69 


0.66 


11.05 


-3.25 


-7.34 


-10.61 


-7.90 


-9.45 


-8.59 


-8.55 


NGC 4742 


4.18 


7.35 


-3.50 


-1.98 


9.80 


-3.05 


-5.35 


-7.34 


-5.48 


-6.85 


-6.11 


-6.10 


NGC 4874 


8.57 


9.64 


-6.07 


1.61 


14.50 


-5.37 


-9.92 


-12.73 


-10.51 


-12.06 


-11.29 


-11.59 


NGC 4889 


8.33 


9.54 


-5.77 


1.60 


14.13 


-5.09 


-9.67 


-12.45 


-10.24 


-11.78 


-11.07 


-11.30 


NGC 524 


6.13 


8.57 


-5.10 


1.00 


12.13 


-4.06 


-8.38 


-11.73 


-8.94 


-10.45 


-9.84 


-9.97 


NGC 5813 


6.44 


8.73 


-5.14 


1.10 


12.51 


-4.29 


-8.67 


-11.95 


-9.21 


-10.79 


-10.09 


-10.21 


NGC 5845 


4.66 


7.69 


-4.41 


0.78 


10.00 


-2.81 


-7.20 


-10.48 


-7.77 


-9.31 


-8.73 


-8.73 


NGC 596 


5.52 


8.24 


-4.61 


0.95 


10.80 


-3.07 


-7.46 


-10.73 


-8.02 


-9.57 


-8.97 


-8.97 


NGC 6166 


8.46 


9.60 


-6.24 


1.69 


14.68 


-5.59 


-10.26 


-13.55 


-10.81 


-12.39 


-11.79 


-11.79 


NGC 720 


6.23 


8.62 


-5.56 


0.96 


13.00 


-4.89 


-9.17 


-11.93 


-9.72 


-11.23 


-10.57 


-10.79 


NGC 7332 


4.70 


7.72 


-3.73 


-0.56 


10.01 


-2.80 


-5.96 


-8.73 


-6.44 


-7.90 


-6.98 


-7.02 


NGC 7768 


7.73 


9.31 


-5.42 


1.62 


13.36 


-4.56 


-9.27 


-12.19 


-9.82 


-11.36 


-10.76 


-10.76 
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Note. — Column 1 is the galaxy name. Column 2 is the V-band mass-to-light ratio and Column 3 gives the corresponding BH mass according to the 
M. — Mb u ige relation llHaring fc Rixll2004h . Columns 4 is the loss cone filling rates for single BHs in non-GR case. Columns 5, 6, and 7 list the stalling 
radius, evolution timescale, and loss cone filling rate for SMBHBs with mass ratio q = 0.1, the results for q = 0.01 and 1 being essentially identical. Both 
t cvo l and ^bi nary are evaluated at a sta n , given by the minimum of equation 1171 The tidal disruption rates for q = 0.01 are given in Columns 8 and 9, 
respectively for the primary and secondary BHs. Columns 10 and 11 have the same meanings as those in Columns 9 and 10, but for q = 0.1, while Columns 
12 and 13 are for q = 1. 
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5.3. Stellar Disruption Rates in the Non-relativistic Case 

First, we ignore GR effect and only consider tidal disruption of solar type stars. In this 
case the stellar disruption rates for single BHs have been given by column 4 in Table [D For 
binaries the rates are calculated according to equation (Q with the loss cone filling rates 
listed in Table [T] and the non-relativistic cross sections from § 14.11 

In Figure [8] we present the stelar disruption rates for both single (triangles) and binary 
(dots and circles) BHs. It is clear that the total disruption rates for SMBHB (N^ inary = 
Nf + Ny) is significantly lower than that for single BH (iV s j le ). The contrast between 
^binary an d ^single i s more than one order of magnitude and increases with q. The significant 
reduction of tidal disruption rates in binary BH systems originates in the depletion of loss- 
cone stars and the suppression of loss cone filling rate by the SMBHB. 

Mechanisms other than two-body relaxation woul d increase the loss cone filling rates 



in both single and binary SMBH systems (MT99; Y02; iMerritt fc Milosavlievic 1120051 ). The 
enhancement in loss cone filling rate is uncertain because it depends on processes which 
are not well understood. However, equation ([7]) implies that the lifetime of hard SMBHB 
decreases with loss cone filling rate. So for long-life SMBHB with th > 10 10 yr the loss cone 
filling rate should not significantly exceed = M 12 /(2J^m,, 10 10 yr). Correspondingly the 
stellar disruption rate in SMBHB system is unlikely to be much greater than 



M d . = F lc . 
cri cri 



Si(r t i) S 2 (r t2 ) 



These upper limits of stellar disruption rates for binary BHs have been presented in Figure M 
as dashed lines. In 48 out of 51 cases iVy nary is below N^ ri , reflecting that two-body relaxation 
is inefficient in refilling the loss cone. For galaxies with M, < 10 8 M Q N^ Ti is always lower 
than N^ in lc , so long-life SMBHBs in these galaxies would manifest themselves by prominently 
suppressing stellar disruption rates no matter what mechanism is responsible for the loss 
cone replenishing. At M. > 1O 8 M could be higher than N!* ngle if M. > 1O 8 M or 
q <C 1, implying that in these galaxies SMBHBs are not so distinct in the efficiency of stellar 
disruption. 



5.4. Stellar Disruption Rates in the Relativistic Case 



In the relativistic case there is a critical mass M„\ above which a BH would swallow 



the whole s tar without tidal disruption (Hills Ill975l ). For solar type stars recent relativistic 
simulation fllvanov fc Chernyakoval 120061 ) gave that M cri ~ (4 — 8) x 1O 7 M if the BH is 
non-spinning and M cri ~ 1O 9 M if the BH is maximally spinning. Due to the ambiguity of 
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Fig. 8. — Tidal disruption rates of solar-type stellar objects as a function of total BH mass 
in Newtonian theory. Dots and circles, respectively, refer to the disruption rates for primary 
and secondary BHs. Triangle are for single BHs. Dashed lines are thresholds above which 
SMBHBs would coalesce within 10 10 yr. The tidal disruption rates are derived under the 
assumptions that each galaxy is spherical and that the loss cone refilling is due to two-body 
relaxation. 
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BH spin it is not clear whether a star approaching a BH more massive than 10 8 M Q would 
necessarily be disrupted and produce a flare. However, it is expected that stellar disruptions 
by these massive BHs are rare because only those stars approaching the spinning BH along 
a corotating orbit close to the equatorial plane could get disrupted. For less massive BHs 
(M. < 1O 8 M ) although equation §1§ suggests that tidal disruption is inevitable (r t > rg), 
part of the stellar debris with low angular momentum ma y directly plunge into the event 



horizon without producing a flare (INolthenus fc Katzlll983l ). The proportion of such debris 



is sensitive to the spin of the star just before tidal disruption and is uncertain so far, but 
the proportion should be low when r ( <C rg. Because of all these uncertainties, we adopted 
a rough assumption that flares are produced only in the condition r t > r$. 

Figure M shows the stellar disruption rates when GR effect is taken into account. The 
tidal disruption rates for binary BHs are calculated with the loss cone filling rates from 
Table [1] and the relativistic tidal disruption cross sections presented in § 14.51 The resulting 
relativistic tidal disruption rates are about r t /(r t — rs) times greater than the non-relativistic 
ones if r t ^ 2rs, or 4rs/r t times greater if rg < r t < 2rg. The dashed lines show the 
thresholds above which SMBHBs would coalesce within 10 10 yr. They are derived 
according to equation (TTgj) but with relativistic cross sections. To study the frequency of 
tidal flares we have set the stellar disruption rates zero when r t < rg, thus the sudden cutoff 
at ~ (1 + g)lO 8 M and ~ + g)10 8 M o is artificial. 

At M. <C 10 8 M Q the stellar disruption rates are similar to those presented in Figured] 
and again we find that the stellar disruption rates for SMBHBs are considerably lower than 
those for single BHs. While at M. ^> 10 s M & the only sources producing tidal flares are the 
secondary BHs, though the rates are low. 



5.5. Effect of a Spectrum of Stellar Masses 

The results of previous sections are based on the assumption that each galaxy is com- 
posed of stars with solar mass and radius. Since none of our sample galaxies exhibits recent 
nuclear star formation, taking into account a mass spectrum of main-sequence stars would 
result in numerous low mass stellar objects which have smaller mean tidal disruption cross 
sections and are less efficient to relax via two-body interaction. But these effects would not 
qualitatively change the tidal disrupti on rates since: (a) the tidal radius is not sensitive to 



the stellar mass because r* ~ mj' 47 (IBond. Arnett fc Carr I Il984l ) so that r t ~ rn® AA ] (b) 



the increment in stellar number density to reproduce the mass distribution inferred from 
observations compensates for the decrease in the efficiency of two-body relaxation. 
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Fig. 9. — Tidal disruption rates of solar-type stellar objects as a function of total BH mass, 
including general relativisitic effects. Symbols and dashed lines have the same meanings as 
those in Figure El 
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However, tidal disruption of off-main-sequence giant stars becomes important when 
BHs are more massive than 10 8 M o because in this case the tidal radius is greater than the 
Schwarzschild radius. Following MT99 we assume that giant stars with time-averaged radius 
r g = 15r contribute g — 1% to the total stellar po pulation. These value s are consistent with 
the approximate stellar evolution model given in ISyer fc Ulmer I Il999i . For single BHs the 
disruption rates of giants in the diffusive loss cone limit are obtained with gF^ ngle (MT99), 
where the loss cone filling rates -F^g^ in relativistic case are from Table [TJ For binary BHs 
the disruption rates are calculated in the same way as in § 15.41 except that the loss cone 
filling rates and tidal radii for giant stars are used. Since we are interested the frequency of 
tidal flares we have set the stellar disruption rates zero when r t < r$. 

Figure [TU] shows the stellar disruption rates and the thresholds N^ ri when both solar 
type and giant stars are taken into account. At M. < 1O 8 M the stellar disruption rates 
are dominated by disruption of solar type stars therefore similar to those in Figure [9j At 
M. > 1O 8 M disruption of giant stars take over in both single and binary BH systems, 
resulting in low but no longer zero flaring rates. The suppression of stellar disruption rate 
by SMBHB is still obvious if the loss cone refilling is dominated by two-body relaxation, but 
other mechanisms could potentially enhance the stellar disruption rates for massive SMBHBs 
(M. > 10 8 M Q ) to a level indistinguishable from those for single BHs. 



6. RATES OF TIDAL FLARES IN LOCAL UNIVERSE 



To estimate the density of tidal flares in local universe we adopted the mass fu nctions 
of SMBH for early (E+SO) and late (Sabcd) type galaxies from iMarconi et al.l (120041 ). which 
are converte d from the distributio n of stellar velocity dispersion using the empirical M. — 
cr* relation (ITremaine et al.l 120021 ) . The mass functions are presented in Figure [TTI and in 
the following we assume that the total BH mass (M12) in SMBHB systems follows the 
same distribution. At M. < 1O 6 M the mass distribution of BHs is unclear because the 
demographics of the BHs in dwarf galaxies is not well established. The total BH mass 
density according to Figure [TTI is 2.5 x 10 5 M Mpc~ 3 (Ho = 70 km s " 1 Mpc -3 throughout 



this paper), consistent with the value given by lYu fc Tremaine I (120021 ). 



The flaring rates per unit volume are calculated by integrating the BH mass functions 
weighted by the flaring rates. The flaring rates for single and binary BHs are from Fig- 
ures (we have smoothed the flaring rates) and the fiducial field of integration is M.(Mi2) G 
[1O 6 M , 1O 1O M ]. In single BH case we do a second set of calculations with M, G [1O 3 M , 10 10 M G 
to account for the flaring rates from the intermediate-mass BHs (IMBHs, 1O 3 M < M. < 
1O 6 M ) at the centers of dwarf galaxies, while for binary BHs we skip such calculations 
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Fig. 10. — Total tidal disruption rates both of stellar-type stellar objects and red giants as a 
function of total BH mass, including general relativistic effect and assuming that red giants 
contribute 1% to the loss cone filling rate. Symbols and dashed lines have the same meanings 
as those in Figure [SJ 
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because most dwarf galaxies have not experienced mergers flHaehnelt fc Kauffmannll2002l ). 
The flaring rates at M m ^ 10 6 M Q are from the linear extrapolation of the logarithmic flaring 
rate at M, > 1O 6 M . We do not consider flares produced by giants because we will see below 
that these flares differ from those produced by normal stars. Our result is not sensitive to 
this treatment because the mass function of SMBHs cuts off steeply above 10 8 M Q where 
flares from giants become dominant. 

The flaring rates per unit volume for different types of galaxies and for both single and 
binary BHs are presented in Table [2J Note that taking into account asymmetric stellar dis- 
tributions and relaxation mechanisms other than two-body interaction would significantly 
increase the loss cone filling rates in both single and binary BH systems, therefore dramat- 
ically increase the flaring rates. According to Table [21 it is obvious that if SMBHBs are 
ubiquitous in the centers of both early and late type galaxies the flaring rate in local uni- 
verse would be more than one order of magnitude lower than that in single BH case. It is 
also clear that dwarf galaxies, if they follow the same mass distribution as that of massive 
galaxies and harbor IMBHs at their centers, would contribute the majority to the flares in 
late type galaxies but have little effect on the flaring rate in early type galaxies. 

If a fraction, /bin, of SMBHs in the investigated galaxies are in binary and a fraction 
1 — /bin are single, the expected tidal disruption rate is 

Nf = (l-f hhl )N[ Qglc + fUN{ + N(). (19) 

N* as a function of /bin is presented in Figure [121 It shows that iV* is not sensitive to q 
unless /bi n — 1. In § [7] we will use these /bin — diagrams to constrain f hin for different 
types of galaxies. 



7. DISCUSSION 

The hierarchical structure formation model in CDM cosmology predicts that SMBHBs 
are continuously formed across the merging history of galaxies. Constraining the abundance 
of SMBHBs among inactive galaxies is essential to test the theoretical models of the dynamic 
evolution of binary BHs. In this paper we have studied the effect of hard binary BHs on 
stellar disruption rates, trying to find distinct observational signatures of SMBHBs in galaxy 
centers. We focus our attention on inactive galaxies which contain the final products of BH 
mergers so we do not consider the effect of gas on SMBHB evolution. 

We have carried out numerical scattering experiments designed for hard SMBHBs to 
investigate the probability of stellar disruption. Since r t <C ah, in each set of experiments 
usually a large number of test particles (N ~ 10 7 ) were used to give statistical meaningful 
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Fig. 12. — Expected tidal disruption rates of stellar objects by black holes as a function of 
the fraction of binaries for all (upper), early (lower left), and late (lower right) type galaxies. 
Solid lines are our calculated flaring rates for spherical galaxies with two-body relaxation. 
Dotted, short-dashed, and dot-dashed lines are calculated by simply multiplying the solid 
lines with, respectively, 2, 10, and 100. The long-dashed horizontal lines indicates the best 
estimates of the flaring rates in real universe and the shaded regions corresponding to an 
assumed statistical uncertainty of one order magnitude. 
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results. For solar type stars the test particle approximation should be valid outside r t as long 
as Eq -C GMi 2 /a. In this situation the self-binding energy of a star ~ Gml/r* is smaller by 
a factor (M./ 'M*) 2 / 3 than its kinetic energy at r t so the loss of orbital energy due to tidal 
dissipation is not important. The test particle approximation may be problematic for giant 
stars because they are much less concentrated and may suffer from mass stripping during 
their encounters with BHs. This uncertainty could be incorporated into the uncertainty in 
r g . In single BH case the stellar disruption rate is not sensitive to r g if the loss cone is 
in the diffusive regime but would be proportional to r g if the loss cone is in the pinhole 
regime (MT99). For binary BHs the loss cone filling rate does not depend on r g but the tidal 
disruption cross section scales with r g when gravitational focusing is important. Therefore 
unless the effective r g deviates significantly from our choice r g = 15r Q our results for giant 
stars would not be qualitatively different. 

We find that in binary BH systems the multi-encounter cross sections are greater than 
the single-encounter ones. This is because binary BHs tend to scatter particles onto tem- 
porary bound orbits which enhances the number of star-binary encounters. If a star is only 
partly disrupted during each close passage from the BHs, multi-flares would be produced 
occur due to the multi-encounters, which provides a distinct observational signature of SMB- 
HBs. However, we have seen in § 14.11 that the probability of multi-flaring may be extremely 
low because r t <C a. Even if r t ~ a equation (fT2]) implies that a/rsi ~ M 8 2 ^ 3 so grav- 
itational wave radiation would quickly drive the BHs to coalesce, making the multi-flares 
unlikely. A secondary flare could also be produced if the secondary member of SMBHB 
captures sufficient debris from the star disrupted by the primary BH, or vise versa. The 
probability of these events depends on the opening angle of all the spewed debris about the 
stellar disrupting BH, which needs to be determined by hydrodynamic simulations. 

We have shown that the normalized cross sections are not sensitive to the star's initial 
binding energy or to the eccentricity of the BH binary, but do depend on the amplitude and 
orientation of the star's initial angular momentum. The cross sections presented in Figures [1] 
and [7J which have been equally averaged over the whole loss cone and over all directions, 
can be directly used if the loss cone of SMBHB is isotropic and is in the pinhole regime. In 
the diffusive loss cone limit the stars falling to the binary BHs satisfy J > J\ c . According to 
Figure H] the differential tidal disruption cross sections at J > J\ c are lower than those inside 
the loss cone, so the mean tidal disruption cross section in the diffusive limit is smaller than 
those presented in Figures [1] and [7J 

We find that the tidal disruption rates in binary BH systems are more than one order of 
magnitude lower than those in single BH systems. The suppression of stellar disruption rate 
by SMBHB originates in the low loss cone filling rate in binary BH systems, which results 
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from the stellar cusp destruction due to the hardening of SMBHBs. Even efficient loss 
cone filling mechanisms are taken into account the contrast between the stellar disruption 
rates in single and binary BH cases are still prominent in less massive galaxies with M. < 
1O 8 M , because the survival of SMBHBs requires that the loss cone filling rate should not 
exceed a critical value which is pro portional to the tot al BH mass of SMBHB. Our result is 
qualitatively different from that of 



Ivanov et al.l (120051 ). who suggested that a SMBHB will 



enhance the stellar disruption rates. The difference is due to that we are studyin g a much 
later evolution stage at which the dense galactic cusps considered by llvanov et al.l have been 
destroyed during the hardening of SMBHBs. According to the calculation in Y02 binary 
BHs reach the hard radius i n about a dynamic f riction timescale of the merging galaxy. So 
the situation consid ered by llvanov et al.l (120051) applies for a short period relative to the 
lifetime of a galaxy. iMerritt fc Wangl (120051 ) suggested that because loss cone refilling takes 
time the stellar disruption rate immediately after SMBHB coalescence is significantly below 
the final steady disruption rate. But this effect is prominent only in massive galaxies with 
M. > 10 8 M Q , where the relaxation timescale is long. Our result implies that in less massive 
galaxies the flaring rate could give stringent constraint to the abundance of SMBHB. 

In larger galaxies with M, > 10 8 M Q tidal flares are dominated by disruption of giants. 
However, their spectral and variable properties may be different from those of solar type 
stars. A tidal flare is expected to be produced at a scale about r t , initially radiating at 
the Eddington limit LecM — 1-3 x 1O 45 (M./1O 7 M ) erg s -1 with a thermal spectrum of 
effective temperature T cff ~ [L Edd / (Anrfa)] 1 ^ = 2.4 x ltf(rJr & )- l / 2 (mjM Q ) l / & Ml 112 K 
and decaying on a timescale t flarc ~ nGM.(2AE)' 3 / 2 ~ 1.1 yr (r ,/r Q ) 3/2 (m, /M ) ~ 1 Ml' 2 , 
where AE ~ GM m r*/r 2 characterizes the span of the binding energy of the debris (IRees 



19881 ). For solar type stars the spectrum of tidal flare peaks at UV or soft X-ray and the 

i li 

luminosity decays in about M 8 year. While for giant stars with r* = 15r and m* = m 

1/2 

the spectrum peaks at optical or UV band and the flare dims on a timescale of about 60M 8 
years. Therefore if a UV/X-ray outburst decaying within months to one year is detected in 
a galaxy with total BH mass much greater than 1O 8 M (determined by other observations), 
it is likely that a less massive secondary SMBH resides in the galactic center. Such event is 
another signature of SMBHB, though Figure [TU1 suggests that the frequency of such event is 
low. 

So far six candidate tidal disruption events have been observed in nearby i nactive galax 
ies. A mong them, one was observed at z = 0.37 by the UV telescope GALEX (IGezari et al. 



20061 ) and the res t were discovere d in galaxies at z < 0.15 during the one year ROSAT All Sky 
Survey (RASS) (IKomossal 120021 ) . The Hubble types for the host galaxies of the five RASS 
flares are not well determined : NGC 5905 is of a SB galaxy an d the other four gala xies look 
like ellipticals (jKomossall2002l ) . From the observations of RASS lDonley et al.l (120021 ) inferred 
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a total flaring rate of N^ hs ~ 2 x 10~ 7 yr^ 1 Mpc -3 for all types of galaxies. From the RASS 
results we can also give a rough estimate of the flaring rates for different types of galaxies: 
^obs ~ ^-^ x y r_1 Mpc -3 for early type galaxies and iV^ bs ~ 7x 10~ 8 yr" 1 Mpc -3 for late 
types. We plot N^ hs as long-dashed horizontal lines in Figure d2] and the corresponding un- 
certainties with shaded areas. Because of the poor statistics in N^ hs and the ill-determination 
of Hubble types of the host galaxies, the results are highly uncertain. 

We did the calculation based on the assumptions that galaxies are spherical and two- 
body scattering relaxation dominates the loss-cone refilling. However, taking into account 
more realistic non-spherical stellar distributions and other loss-cone refilling mechanisms 
besides two-body interaction will probably significantly increase the loss-cone refilling rates 
in both single and binary SMBH systems. MT99 and Y02 showed that taking into account 
axisymmetric stellar distribution will increase the loss cone refilling rates by a factor of 2 for 
both single SMBHs and hard SMBHBs. The loss cone refilling rates in triaxial galaxies are 
slightly higher than those in axisymmetric ones if the triaxiality is les s than 0.1 but can be tw o 
orders of magnitude higher if the triaxiality is more extreme (Y02; iMerritt fc Poon 1 120041 ) . 
Resonant relaxation would increase the loss cone refilling rates by a factor of 2 — 10 relative 
to those in the case of two-body relaxation, depending on the eccentriciti es of stellar orbits 
( IRauch fc Tremainelll996l ; iRauch fc Ingallslll998l ; iGiirkan fc Hopmanl|2007l ). To calculate the 
tidal disruption rates by taking into account the above physical processes is out of the scope 
of this paper. To illustrate the effects of increased loss-cone refilling rates and tidal disruption 
rates on the estimated fraction of binary BHs, /bin, we simply multiply our calculated stellar 
disruption rates N* by 2, 10, and 100 and plot the results in Figure [T2J In Figure [T21 the 
intersections of long-dashed and solid lines give the lower limit to /bin, while the intersections 
of long-dashed and dot-dashed lines roughly give the upper limits. When all types of galaxies 
are considered, binary fraction / bin is between 0.4 and 1.0. For early type galaxies, the solid 
lines do not intersect the long-dashed line, implying that other loss-cone filling mechanisms in 
addition to two-body relaxation should also be important. For late type galaxies, Figure [T21 
suggests a binary fraction /bi n ~ 75% and that extreme triaxiality may not be common, 
otherwise the current observation would imply an extremely high fraction of SMBHBs. We 
would like to note that the constraints on /bin are very uncertain due to the large uncertainty 
in the current iV^ bs . Detection of more candidate tidal disruption events is needed to give 
better constraints on /bin and on the dominant loss-cone filling mechanisms in different types 
of galaxies. 



The non-detection of X-ray or UV flares in dwarf galaxies seems inconsistent with the 
calculated high flaring rates. The discrepancy may be due to (1) a lack of sufficient sensitivity 
in current surveys, (2) the continuous a nd steady accretion of st ellar debris onto IMBHs 
during successive tidal disruption events (jMilosavljevic et al.ll2006l ). or (3) a lack of IMBHs 
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in the centers of dwarf galaxies. Future X-ray and UV surveys with improved sensitivity and 
angular resolution would help constraining the fraction of IMBHs in galaxy centers and give 



interesting implicatio ns to the formation and the merging history of IMBHs fjVolonteri et al. 



200alVolonterill2007f ). 



8. CONCLUSIONS 

To investigate whether SMBHBs are ubiquitous in nearby inactive galaxies or coalesce 
rapidly in galaxy mergers, we have studied the tidal disruption rates of stellar objects in 
both single and binary SMBH systems. We have calculated the interaction cross sections 
between hard SMBHBs and intruding stars by carrying out intensive numerical scattering 
experiments with typically 10 7 particles, taking into account the initial binding energy and 
angular momentum of particle, the eccentricity of the orbit of SMBHB, and including general 
relativistic effects. We have also derived empirical formulae for the relativistic cross sections, 
which can be applied to SMBHBs with a wide range of semimajor axis and todal BH mass. 

We have calculated the rate of loss cone refilling due to two-body relaxation for a 
sample of 51 nearby galaxies, assuming that each galaxy is spherical. The steady loss cone 
filling rates in binary BH systems would be significantly suppressed due to the three-body 
interaction between SMBHBs and stars passing by. We have calculated the tidal disruption 
rates respectively for single and binary SMBHs by combining the loss cone filling rates and 
the tidal disruption cross sections. We find that the tidal disruption rate in SMBHB systems 
is more than one order of magnitude lower than that in single SMBH systems. For galaxies 
with BHs more massive than 1O 8 M , a UV/X-ray flare at galactic center decaying within 
one year provides strong evidence of a secondary BH, although the probability of such events 
is low. 

Finally we have calculated the flaring rates in local universe using the BH mass function 
given in the literature. The comparison of the calculated flaring rates and the preliminary 
results from current X-ray surveys could not yet tell whether SMBHBs are ubiquitous in local 
universe. Future UV/X-ray surveys with improved sensitivity and duration are needed. 

We are grateful to Dr. Xue-Bing Wu, Dr. Stefanie Komossa, and Dr. Vladimir Karas for 
many constructive comments. We thank the referee for helpful comments and suggestions. 
Many thanks are due to Bing-Xiao Xu, Ran Wang, Lei Qian and Zhao-Yu Li for fruitful 
discussions. During this work we have used the SGI Altix 330 system at the Department 
of Astronomy, Peking University (PKU) and CCSE-I HP Cluster of PKU. This work is 
supported by the National Natural Science Foundation of China (No. 10573001) and by 



-39- 



the national 973 program (No. 2007CB815405). JM thanks the Royal Society for financial 
support. 

REFERENCES 

Armitage, P. J. & Natarajan, P., 2002, ApJ, 567, L9 

Begelman, M. C, Blandford, R. D., & Rees, M. J., 1980, Nature, 287, 307 
Berczik, P., Merritt, D., Spurzem, R., & Bischof, H.-P, 2006, ApJ, 642, L21 
Bond, J. R., Arnett, W. D., & Carr, B. J., 1984, ApJ, 280, 825 
Chatterjee, P., Hernquist, L., Loeb, A., 2003, ApJ, 592, 32 
Cohn, H. & Kulsrud, R. M., 1978, ApJ, 226, 1087 
Croton, D. J. et al, 2006, MNRAS, 365, 11 

Donley, J. L., Brandt, W. N., Eracleous, M., & Boiler, T., 2002, AJ, 124, 1308 

Escala, A., Larson, R. B., Coppi, P. S., & Mardones, D., 2005, ApJ, 630, 152 

Ferrarese, L. & Ford, H., 2005, SSRv, 116, 523 

Ferrarese, L. & Merritt, D., 2000, ApJ, 539, L9 

Gaskell, C. M., 1985, Nature, 315, 386 

Gebhardt, K. et al., 2000, ApJ, 539, L13 

Gezari, S. et al., 2006, ApJ, 653, L25 

Gould, A. & Rix, H. W., 2000, ApJ, 532, L29 

Giirkan, M. A. & Hopman C, 2007, MNRAS, 379, 1083 

Haehnelt, M. G. & Kauffmann, G., 2002, ApJ, 336, 61L 

Hairer, E., Norsett, S. P., & Wanner, G., 1987, Solving Ordinary Differential Equations I 
(1st ed.; Berlin: Springer- Verlag) 

Halpern, J. P., Gezari, S. & Komossa, S., 2004, ApJ, 604, 572 

Haring, N. & Rix, H.-W., 2004, ApJ, 604, 89L 



-40- 

Hernquist, L. & Mihos, J. C, 1995, ApJ, 448, 41 
Hills, J. G., 1975, ApJ, 254, 295 
Hills, J. G., 1983, ApJ, 88, 1857 

Hopkins, P. F., Somerville, R. S., Hernquist, L., Cox, T. J., Robertson, B., & Li, Y., 2006, 
ApJ, 652, 864 

Hut, P., 1983, ApJ, 272, L29 

Ivanov, P. B. & Chernyakova, M. A., 2006, A&A, 448, 843 

Ivanov, P. B., Papaloizou, J. C. B., & Polnrev, A. G., 1999, MNRAS, 307, 79 

Ivanov, P. B., Polnarev, A. G., & Saha, P., 2005, MNRAS, 358, 1361 

Kauffmann, G. & Haehnelt, M. G., 2000, MNRAS, 311, 576 

Kazantzidis, S. et al., 2005, ApJ, 623, L67 

Komossa, S., 2002, Rev. Mod. Astron., 15, 27 ( |astro-pli/0209007| ) 

Komossa, S., 2006, MmSAI, 77, 733 

Komossa, S. & Greiner, J., 1999, A&A, 349, L45 

Komossa, S., Halpern, J., Schartel, N., Hasinger, G., Santos-Lleo, M., & Predehl, P., 2004, 
ApJ, 603, L17 

Liu F. K., 2004, MNRAS, 347, 1357 

Liu, F. K. & Chen, X., 2007, ApJ, 671, no. 2 (astro-ph/0705.1077) 

Liu, F. K., Liu, B. F., & Xie, G. Z., 1997, A&AS, 123, 569 

Liu, F. K. & Wu, X. B., 2002, A&A, 388, L48 

Liu, F. K., Wu, X. B., & Cao, S. L., 2003, MNRAS, 340, 411 

Liu, F. K., Xie, G. Z., & Bai, J. M., 1995, A&A, 295, 1 

Liu, F. K., Zhao, G., & Wu, X. B., 2006, ApJ, 650, 749 

Lu, Y., Yu, Q., & Lin, D. N. C, 2007, ApJ, 666, L89 



-41 - 

Marconi, A., Risaliti, G., Gilli, R., Hunt, L. K., Maiolino, R., & Salvati, M., 2004, MNRAS, 
351, 169 

Magorrian, J. et al, 1998, AJ, 115, 2285 

Magorrian, J. & Tremaine, S., 1999, MNRAS, 309, 447 (MT99) 

McLure, R. J. & Dunlp, J. S., 2002, MNRAS, 331, 795 

Merritt, D., 2006, ApJ, 648, 976 

Merritt, D. & Ekers, R. D., 2002, Sci, 297, 1310 

Merritt, D. & Milosavljevic, M., 2005, Liv. Rev. Rel, 8, 8 

Merritt, D. & Poon, M. Y., 2004, ApJ, 606, 788 

Merritt, D. & Wang, J.-X., 2005, ApJ, 621, L101 

Milosavljevic, M. & Merritt, D., 2001, ApJ, 563, 34 

Milosavljevic, M. & Merritt, D., 2003, ApJ, 596, 860 

Milosavljevic, M., Merritt, D. & Ho, L. C, 2006, ApJ, 652, 120 

Nolthenus, R. A. & Katz, J. I., 1983, ApJ, 269, 297 

Ostorero, L., Villata, M., & Raiteri, C. M., 2004, A&A, 419, 913 

Paczyriski, B. & Witta, P. J., 1980, A&A, 88,23 

Peters, P. C, 1964, Phys. Rev. B, 136, 1224 

Quinlan, G. D., 1996, New A, 1, 35, (Q96) 

Raiteri, C. M. et al, 2001, A&A, 377, 396 

Rauch, K. P. & Ingalls, B., 1998, MNRAS, 299, 1231 

Rauch, K. P. & Tremaine, S., 1996, New A, 1, 149 

Rees, M. J., 1988, Nature, 333, 523 

Richstone, D. O. et al. 1998, Nature, 395, 14 

Romero, G. E., Chajet, L., Abraham, Z. & Fan, J. H., 2000, A&A, 360, 57 



-42- 



Sillanpaa, A., Haarala, S., Valtonen, M. J., Sundelius, B., & Byrd, G. G., 1988, ApJ, 325, 
628 

Sesana, A., Haardt, F., & Madau, P., 2006, ApJ, 651, 392 
Springel, V. et al., 2005, Nature, 435, 629 
Syer, D. & Ulmer, A., 1999, MNRAS, 306, 35 
Tremaine, S. et al, 2002, ApJ, 574, 740 
Ulmer, A., 1999, ApJ, 514, 180 

Valtaoja, E., Terasranta, H., Tornikoshi, M., Sillanpaa, A., Aller, M. F., Aller, H. D. & 
Hughes, P. A., 2000, ApJ, 531, 744 

Villata, M. & Raiteri, C. M., 1999, A&A, 347, 30 

Volonteri, V., 2007, ApJ, 663, LI 

Volonteri, M., Haardt, F., & Madau, P., 2003, ApJ, 582, 559 

Wang, J.-X. & Merritt, D., 2004, ApJ, 600, 149 

Yu, Q., 2002, MNRAS, 331, 935 (Y02) 

Yu, Q. & Tremaine, S., 2002, MNRAS, 335, 965 

Zier, G, 2005, MNRAS, 364, 583 



This preprint was prepared with the AAS IATgX macros v5.2. 



-43- 



Table 2. Flaring rate per unit volume in yr 1 Mpc 



Galaxy type 


N f , 

single 




single 




Q 




n{ 






n( 




(1) 


(2) 




(3) 




(4) 




(5) 






(6) 




Early 


6.0 x 10" 


-8 


6.0 x 10' 


-8 


0.01 


3.6 


x 10- 


-10 


3.8 


x 10- 


-13 


Early 


6.0 x 10" 


-8 


6.0 x 10" 


-8 


0.1 


1.2 


x 10- 


-10 


2.4 


x 10- 


-12 


Early 


6.0 x 10" 


-8 


6.0 x 10- 


-8 


1 


2.1 


x 10- 


-11 


2.1 


x 10- 


-11 


Late 


3.0 x 10' 


-7 


1.4 x 10- 


-6 


0.01 


3.8 


x 10' 


-9 


2.2 


x 10- 


-11 


Late 


3.0 x 10' 


-7 


1.4 x 10- 


-6 


0.1 


2.3 


x 10' 


-9 


9.0 


x 10- 


-11 


Late 


3.0 x 10' 


-7 


1.4 x 10- 


-6 


1 


5.6 


x 10' 


-10 


5.6 


x 10- 


-10 



Note. — Column 1 gives the galaxy type: 'Early'=E+S0, 'Late'=Sabcd. 
Columns 2 and 3 are the flaring rates per unit volume in single BH case, the 
former being calculated for M, > 1O 6 M while the latter for M. > 1O 3 M . 
q is the assumed mass ratio of SMBHB, and N( and are respectively 
the flaring rates for primary and secondary BHs in binary BH case. 



